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ABSTRACT 

We have investigated the broad-hne Balmer decrements (Ha/H/3) for a large, homo- 
geneous sample of Seyfcrt 1 galaxies and QSOs using spectroscopic data obtained in 
the Sloan Digital Sky Survey. The sample, drawn from the Fourth Data Release, com- 
prises 446 low redshift {z < 0.35) active galactic nuclei (AGN) that have blue optical 
continua as indicated by the spectral slopes in order to minimize the effect of dust ex- 
tinction. We find that (i) the distribution of the intrinsic broad-line Ha/H/3 ratio can 
be well described by log-Gaussian, with a peak at Ha/H/3 =3.06 and a standard devi- 
ation of about 0.03 dex only; (ii) the Balmer decrement does not correlate with AGN 
properties such as luminosity, accretion rate, and continuum slope, etc.; (iii) on aver- 
age, the Balmer decrements are found to be only slightly larger in radio-loud sources 
(3.37) and sources having double-peaked emission-line profiles (3.27) compared to the 
rest of the sample. We therefore suggest that the broad-linc Ha/H/3 ratio can be used 
as a good indicator for dust extinction in the AGN broad-line region; this is especially 
true for radio-quiet AGN with regular emission-line profiles, which constitute the vast 
majority of the AGN population. 
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1 INTRODUCTION 

Hydrogen Balmer decrements are often used to determine 
the amount of dust extinction attenuating the observed 
emission lines because the intrinsic decrements of Balmer 
recombination lines are quite insensitive to the gas temper- 
ature and density in low density, dilute radiation field con- 
ditions (Osterbrock 1989). Particularly, the Ha/H/? ratio is 
most frequently used for the strongness of the Ha and H/3 
lines and their relatively large wavelength span, which ren- 
der the derived amount of extinction less affected by mea- 
surement uncertainties. For HH region photoionized by a hot 
star, an intrinsic Ha/H/3 value of 2.87 is found, as predicted 
by Case B recombination (at a typical electron density of 
< 10* cm~^ and temperature of 10* K). A value of 3.1 is 
generally adopted for the narrow line region (NLR) of ac- 
tive galactic nuclei (AGN), where Ha emission is slightly 
enhanced by coUisional excitation due to the presence of 
gas of higher densities and the presence of partly ionized 



transition region resulting from much harder ionizing con- 
tinuum (Gaskell & Ferland 1984, Halpern & Steiner 1983). 
For the broad line region (BLR) of AGN, however, where 
the density is so high (typically > 10^ cm"'') that col- 
lisional, optical-depth and radiative-transfer effects become 
important as predicted by BLR photoionization modeling 
(e.g., Netzer 1975, Kwan & Krolik 1981, Canfield & Puetter 
1981, CoUin-Souffrin et al. 1982, Rees et al. 1989, Dumont 
et al. 1998, Korista & Goad 2004, see Baldwin 1997 for a re- 
view I , the broad-line Ha/H/? ratio (HaVH/3'') varies widely 
in different BLR conditions. Observationally, the Ha^'/H/?'' 
ratios are larger (steeper) than the Case B recombination 
value in most Seyfert Is and QSOs (e.g., Osterbrock 1977, 
Wu et al. 1980, Rafanelli 1985, and Fig. 3 of Dong et al. 
2005); moreover, in broad-line radio galazxies and Seyfert 
1.8/1.9 galaxies HaVH/?'' can be as steep as 10 or higher 
(Osterbrock et al. 1976, Osterbrock 1981, Crenshaw et al. 
1988). Thus, although it is unclear yet whether the observed 
steep Ha''/H/3'' is due to the high-density effects mentioned 
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above, or just to wavelength-dependent extinction by dust 
(e.g., Osterbrock 1984, Goodrich 1995), it has been gener- 
ally believed that the Ha^'/H/?'' ratio cannot be used as an 
indicator of reddening in BLR, as is in AGN NLRs or HII 
regions. 

However, we have noticed that the range of the 
Hq^'/H/?'' ratios of Seyfert 1/QSOs is fairly small, typically 
from 2.5 to 5 (see references above), probably suggestive 
of a small dispersion of the intrinsic Hq''/H/3'' ratio. For in- 
stance, for 94 Seyfert Is and QSOs having u'—g' < 0.6 culled 
from the Sloan Digital Sky Survey (SDSS; York et al. 2000) 
Early Data Release (EDR; stoughton et al. 2002), the stan- 
dard deviation of Ha/H/3 is only 0.36 around a mean of 2.98 
(Dong et al. 2005). This fact motivated us to explore sys- 
tematically the intrinsic Balmer decrements in AGN, their 
distribution and potential dependence on other AGN prop- 
erties (e.g. radio-loudness, luminosity, accretion rate, line- 
profile parameters, continuum slope), by taking advantage of 
the unprecedented spectroscopic data from the SDSS. This 
study will be able to address the question as to whether 
the HqVH/3'' ratio can be taken as an indicator of the BLR 
extinction, at least in a statistical manner and for specific 
sub-classes of AGN. To this end, we have compiled a larger 
sample of blue AGN, including both Seyfert Is and QSOs, 
from the SDSS Fourth Data Release (Adelman-McCarthy et 
al. 2006) as to the shape of their continua. We assume a cos- 
mology with Ho=70 km s"^ Mpc"\ and riA=0.7. 



2 SAMPLE CONSTRUCTION AND DATA 
ANALYSIS 

2.1 Sample Definition 

Our aim is to select blue AGN that are free of dust ex- 
tinction. It has been noted that even in quasars with rel- 
atively blue colors, such as Palomar-Green (PG) quasars 
{U-B < -0.44, Schmidt & Green 1983), there is noticeable 
internal dust extinction inside the AGN (Rowan- Robinson 
1995, Baskin & Laor 2004). A color criterion bluer than 
the average of QSOs is desirable. The average optical-near- 
ultraviolet slope of QSO continuum is found to be a\ ~ 1.5 
{fx = A~"-^ ; Vanden Berk et al. 2001 and references listed 
in their Table 5). We thus define blue AGN as those with 
a continuum slope a\ > 1.5, where a\ is fitted in the rest- 
wavelength range of 4030-5600 A. We consider thus selected 
AGN to be least affected by dust extinction. This point is 
to be further discussed in §4.1, as well as the representative- 
ness of our sample. In practice, we limit redshifts z < 0.35 so 
that the Ha line lies within the SDSS spectral coverage. We 
select objects with the median spectral signal-to-noise ratio 
(S/N) ^ 10 per pixel only to ensure accurate measurement 
of the Balmer decrement. 



2.2 Overview of Data Processing 

Applying the above redshift and S/N cutoffs yielded a pool 
of ~4100 objects classified as AGN in the SDSS DR4 spec- 
tral data set, from which our sample is to be culled. The 
spectra are first corrected for Galactic extinction using the 
extinction map of Schlegel et al. (1998) and the reddening 
curve of Fitzpatrick (1999). The spectra are transformed 



into the rest frame using the redshift as determined from the 
peak of the [O IIIJA5007 emission line. In order to measure 
accurately broad Balmer lines, we have to subtract prop- 
erly the continuum, the Fell emission multiplets, and other 
emission lines nearby. In common practice, the subtraction 
is performed step by step: first to fit and subtract the AGN 
continuum and then Fell multiplets (or the opposite) and 
finally to fit emission lines (e.g., Boroson & Green 1992, 
Marziani et al. 2003); or first to fit and subtract simultane- 
ously the continuum and the Fell emission multiplets (so- 
called "pseudo-continuum") and then to fit other emission 
lines (e.g., Dong et al. 2005, Greene & Ho 2005, Zhou et al. 
2006). Unfortunately, for the optical spectra of most Seyfert 
Is and QSOs, fitting the continuum and the Fell emission is 
highly complicated by several facts as follows. 1) There are 
essentially no emission-line-free regions where the contin- 
uum can be determined (e.g., Vanden Berk et al. 2001). 2) 
The Fell AA4434 — 4684 features, generally prominent, are 
often blended with broad lines of H7, Hell A4686 and H/3. 3) 
Often the QSO continuum cannot be described by a single 
power-law from H(5 to Ha, which means that we have to de- 
termine the local continuum for the H/3 and Ha regions sepa- 
rately. Limited by these complications, the common step-by- 
step spectral subtraction procedure cannot achieve Balmer 
decrement measurement accurate enough for our purpose 
in this work, which is rather sensitive to the measurement 
uncertainties of the emission-line fluxefl Here we adopt an 
alternative method to fit simultaneously the continuum, the 
Fe II and the other emission lines, giving emphasis on proper 
determination of the local pseudo-continua. If there are still 
large residuals left in the H/3 and/or Ha regions, a refined fit 
of the emission-lines is performed to the pseudo-continuum 
subtracted spectra. 

Having completed the spectral fitting procedure, we se- 
lect blue AGN according to the above slope criterion. Several 
objects having many bad pixels in the H/3 or Ha regions are 
removed. For eleven objects having duplicated spectra, we 
retain the one tagged as 'SciencePrimary' only according 
to the SDSS Catalog Archive ServeiQ. Our final sample is 
composed of 446 Seyfert Is and QSOs. Our spectral fitting 
method is described in detail in §2.3 and 2.4. 



2.3 Simultaneous Fit of Continuum and Emission 
Lines 

As discussed above, we need to fit simultaneously the re- 
spective local continua in the H/3 and Ha regions and the 
Fe II emission spectrum. We follow the procedure described 
in below, which is implemented using IDL routines[f| 

1. We fit each SDSS spectrum in the rest- wavelength 
range from 4030A to 7500A, assuming a broken power law 
with a break wavelength of 5600A, i.e., aiA~"^-^ for the HtJ- 
H7-H/3 region and a2X °''^-^ for the Ha region. We limit the 



^ For instance, with a typical error of 10% for the omission-line 
fluxes as in Osterbrock 1977, the error of the Ha/H/3 ratios is 
14% (0.061 dex), which is fairly large for our purpose (cf. §3.1). 
^ http://cas.sdss.org/ 

We use the MPFIT package for nonlinear fitting. MP- 
FIT is kindly provided by Craig B. Markwardt, available at 
|http://cow.physics.wisc.edu/~craigm/idl7] 
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wavelength range redward of 4030A to avoid the broad-hne 
He emission. 

2. The optical Fell emission is modeled as C(A) = 
C5C(,(A) + c„C„(A), where Cb{\) represents the broad Fell 
lines with the relative intensities fixed at those of I ZW 1 as 
given in Table A.l of Veron-Cetty et al. (2004), and with 
the same line profile as the broad H/3 line. The redshift of 
the broad Fell lines relative to H/3 is fitted as a free pa- 
rameter. Cn(A) denotes the narrow permitted and forbidden 
Fe II lines with their relative intensities fixed at those listed 
in Table A.2 of Veron-Cetty et al (2004); which have the 
same redshifts and line profiles as the narrow H/3 line. Q 

3. Emission lines other than Iron lines identified from 
the composite SDSS QSO spectrum (see Table 2 in Vanden 
Berk et al. 2001) from H5 to [S II] AGTSlE] are modeled as 
follows. Broad Hydrogen Balmer lines are assumed to have 
the same redshifts and profiles, and each is modeled with 
1-4 Gaussians. The broad Hell A4686 line is modeled with 
one Gaussian. [O HI] A4363 and the AA4959, 5007 doublet 
are assumed to have the same redshifts and profiles, and 
each is modeled with 1-2 Gaussians. Other narrow lines are 
modeled with one Gaussian. Narrow Balmer lines, [N II] and 
[S II] doublets are assumed to have the same redshift and 
profile. The fiux ratio of the [Nil] doublet A6583/A6548 is 
fixed to the theoretical value of 2.96; the [O HI] doublet is 
similarly constrained. 

We note that in many cases the fit is not good in several 
wavelength regions, such as the small region around the min- 
imum between Fell 37,38 complex and H/3, the joint part 
of H/3 and [O III] A4959. These are likely caused by the im- 
perfection of the models, for instance, over-estimation of the 
continuum between the Fell 37,38 complex and H/3, over- 
fit of H/3 with a spurious very broad component, over-fit of 
[OIII] with a spurious blue-shifted wing. These problems 
can be solved by assigning additionally larger weights to 
these critical regions in the fitting, following the weighting 
methods adopted by Tran et al. (1992) and Reichard et al. 
(2003). The exact weights are determined by trial-and-error 
as those which give the best fits — with the minimum cal- 
culated using the original weights (errors) — in the relevant 
regions. 

We also note that the broad lines of He I A4922 and 
A5016 may contribute to the "red shelf" of H/3 (e.g., Veron et 
al. 2002). To investigate this possibility, we compare the H/3 
profiles with the Ha profiles in our sample. Only 4 per cent 
of the sources are found to have apparent redward excess of 
H/3, and the strength is typically less than 5 per cent of H/3. 
We thus consider the contribution of the potential He I lines 
to be negligible. 



Instead of the template spectrum of the broad Fe II emission of 
IZWl readily provided by Veron-Cetty et al. (2004), we use two 
sets of Fell emission templates in analytical forms constructed 
from their measurements, one for the broad line system LI and 
the other for low-excitation narrow line system N3. In addition, 
we also add into our templates Till, Nill, CrII lines listed in 
their Table A.l and A.2. 

^ Several emission-line regions in the 4030-7500A range arc sim- 
ply masked out, because either they are too weak to constrain in 
the fit or they have little effect on the results. These include He I 
A4471, [NI] A5200, [Ca V] A5310, [CI III] A5538, He I A5876, Hcl 
A7066, [Ar III] A7136 and [O II] A7320. 



The assumption that the broad Balmer lines have the 
same profiles is useful to well constrain these lines in the 
fitting since they are highly blended with other lines nearby. 
This assumption is found to have little effect on the deter- 
mination of the local continua (see §2.2), although it is well 
known that broad Balmer lines have slightly different pro- 
files (e.g. Osterbrock & Shuder 1982). A similar situation 
arises in the comparison of the [OIIIJA4363 and A5007 line 
profiles, as the former being often broader. 

We notice that in many objects the relative intensities 
of various Fe II multiplets, such as 37*33 , are more or less 
different from those of IZWl, resulting apparent residuals 
left in some Fell multiplet regions. For about two dozen 
such spectra we re-fit them with different scaling factors 
for the Fell emission blueward and redward of H/3. This 
yields better matched Fe II emission spectrum; the local con- 
tinua and the H/3 and Ha fiuxes remain unchanged, because 
the wavelength regions determining the continua and the 
"H5-f H7-(-H/3-f Ha" lines have much larger total weights 
than that determining the abnormal Fell multiplets. We 
thus adopt the fits that used the uniform factor for the whole 
Fe II spectrum, since we are not concerned with the proper- 
ties of Fell emission in this work. Example spectra as well 
as the fitting results are demonstrated in Fig. 1. 

2.4 Refined Fitting of Emission-line Profile 

We re-calculate the reduced of the fits around H/3 and 
the Ha regions (4750-5050A and 6400-6800A, respectively) 
using the original errors for the 446 spectra. Spectra having 
a fit with the reduced > 1-1 are picked up for further 
refined line-profile fitting using the code described in detail 
in Dong et al. (2005) . Briefiy, we fit each pseudo-continuum 
subtracted spectrum using various schemes and the one with 
the minimum reduced x^ is adopted. Broad lines are fitted 
with multiple Gaussians, as many as 4 at most; the fits are 
accepted when the reduced x^ ^5 1-1 or it cannot be im- 
proved significantly by adding in one more Gaussian (up to 
4) with a chance probability less than 0.05 according to F- 
testlfl Narrow Balmer lines mostly have similar profiles to 
[Nil], [SII] or the line core of [OIII] A5007. At this stage, 
the broad Balmer lines are not required to have the same 
profiles. Fig. 2 shows examples of refined line-profile fitting 
in the H/3 and Ha regions. The line fluxes of H/3 and Ha 
are listed in Table 1, and the broad- line profile parameters 
are listed in Table 2. The data and the fitted spectral pa- 
rameters are available online for the decomposed spectral 
components (continuum. Fell and other emission lines) for 
the 446 objectsQ 

Our two-step procedure of pseudo-continuum and line- 
profile fitting has proved to be robust and self-consistent. 
As a reliability check, we perform the first-step fitting (as in 
§2.3) by setting the initial values of the free parameters of 

^ We have found through our experiment that these criteria 
based on x'^-test and F-test work well, although, theoretically, 
these goodness-of-fit tests holds only for linear models (cf. Lup- 
ton 1993). 
^ Available at 

http: / /staf f. ustc.edu.cn/~xbdong/Data_Release/blueAGN_DR4/[ 
together with an auxiliary code to explain the parameters and to 
demonstrate the fitting. 
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emission lines (except Iron lines) to those yielded from the 
refined line-profile fitting, and find that the best-fits of both 
the continuum and the Fell emission are almost unchanged. 



2.5 Estimation of Parameter Uncertainties 

As we have argued above, estimation of the measurement 
errors of emission-line fiuxes is important for deriving the in- 
trinsic distribution of the Balmer decrements. The errors of 
line fluxes provided by MPFIT are unreasonably small: the 
median errors of broad H(5, H7, H/3 and Ha are 13%, 12%, 
6% and 4%, respectively. They do not account for the uncer- 
tainty introduced by the pseudo-continuum subtraction that 
severely complicates measurement of the broad Balmer lines 
(see Marziani et al. 2003 for a detailed discussion). To take 
this and other possible effects into account, we adopt a boot- 
strap approach to estimate the typical errors for the whole 
sample. We generate 500 spectra by randomly combining the 
scaled model emission lines of one object (denoted as 'A') 
to the emission-line subtracted spectrum of another object 
(denoted as 'B'). The emission line model of object 'A' is 
scaled in such a way that it has the same broad H/3 flux as 
object 'B', in order to minimize changes in S/N within the 
emission-line spectral regions in the simulated spectra. Then 
we fit the simulated spectra following the same procedure as 
described in §2.3 and §2.4. For each parameter, we consider 
the error typical of our sample to be the standard deviation 
of the relative difference between the input {xi) and the re- 
covered {xo), , These relative differences turn out to 
be normally distributed for every parameter. The thus es- 
timated typical Icr errors for the fluxes of broad H5, H7, 
H/3 and Hq are 19%, 19%, 8% and 5%, respectively. Hence 
the measurement error of the Ha''/H/3'' ratio follows the log- 
normal distribution with lcr=0.04 dex. The uncertainties for 
the broad H(5 and H7 line fluxes are large because they are 
relatively weak and/or are often severely contaminated by 
the Fell and [OHIJA4363 emission. Considering the large 
uncertainties of their measurements, we do not discuss the 
H5VH/3'' and H7VH/3'' ratios in this study. 

The errors provided by MPFIT are unreasonably small 
also for the power-law indices and normalizations, commonly 
being only 0.5 and 0.4 per cent, respectively. Using the above 
boot-strap approach, the typical errors are 8 per cent for 
qa,! and Qa,2, 5 and 3 per cent for ci and C2, respectively. 
The mean errors of the fluxes of narrow H/3 and Ha are 6 
and 5 per cent, respectively, which are almost unaffected by 
the pseudo-continuum subtraction. 



3 RESULTS 

3.1 Intrinsic broad line Ha/H/3 Distribution 

The measured Ha''/H/3'' ratios range from 2.3 to 4.2, with 
a mean of 3.1 and skewed towards the large ratio end. Con- 
sidering the log-normal distribution of their measurement 
errors, we plot the histogram in base-10 logarithm form by 
dividing the sample into 20 bins, as shown in Fig. 3. The 
profile of the distribution is very similar to Gaussian. We 
fit the distribution with a Gaussian function by minimiz- 
ing assuming Poissonian errors for the counts in each 



bin. A good fit (see Fig. 3) is achieved with a minimum re- 
duced = 0.89 (17 de grees of freedom) , yielding a mean of 
0.486±0.002 and a standard deviation of 0.046±0.002. This 
model distribution is identical to the histogram at the 96.5% 
significance level according to the Kolmogorov-Smirnov test. 
We re-fit the histogram by varying the number of bins from 
10 through 30 and find that all the fits yield similar re- 
sults at similar confidence levels. Therefore, the observed 
distribution of the Hq''/H/3'' ratios is well described by log- 
normal, with a peak at Hoi'/Hfi^ = 3.06. This fact, to- 
gether with that the measurement errors have a log-normal 
distribution (see §2.5), indicate that the distribution of the 
intrinsic Ha'/H/?'' ratios should be, at least very close to, 
log-normal. The intrinsic Hq'/H/?'' distribution can be ap- 
proximated by deconvolving the measurement uncertainty 
(Icr ~ 0.04 dex) from the observed distribution, yielding a 
standard deviation as small as about 0.03 dex. The decon- 
volved intrinsic distribution is also displayed in Fig. 3. 

It should be noted that measurement of the broad Ha 
and H/3 line fiuxes is insensitive to the deblending with the 
narrow Ha and H/3 components because of the weakness 
of the latter in blue AGN, which account for only ~5 and 
~4 per cent of the total Ha and H/3 fluxes in the sample. 
The distribution of the total line flux ratios log( ^^totai ) hs^s 
a mean of 0.493 and a standard deviation of 0.046, nearly 
identical to those of the broad components log(-^2j.). 



3.2 Dependence of HaVH/3'' on AGN Properties 

We conduct comprehensive statistical analysis using the 
present sample in an attempt to investigate whether there 
exist any systematic trends (biases) of the Ha''/H/3'' de- 
pendence on AGN types and/or properties, such as radio- 
loud objects, objects with clumpy /double-peaked line pro- 
files, objects in low states. This is important as any bi- 
ases, if being significant, would affect the observed distribu- 
tion of Ha''/H/3''. Correlation analysis is performed between 
Ha''/H/3'' and various AGN properties, and the results are 
summarized in Table 3. We report the Spearman rank cor- 
relation coefficient (p) and the probability (Pnuii) that a cor- 
relation is not present. When upper limits are present, we 
use the generalized Spearman rank correlation test as imple- 
mented in the ASURV package (Isobe et al. 1986). We also 
perform statistical tests to compare the Balmer decrements 
among various subsamples. 



3.2.1 Luminosity and Eddington Ratio 

We calculate monochromatic continuum luminosity I/5100 = 
\L\ at 5IOOA from the power-law fit a\\~°'^-^ . Due to the 
blue slope criterion the optical luminosities have little con- 
tribution from host galaxy starlight. Based on the empiri- 
cal BLR radius-luminosity (R-L) relationship and the as- 
sumption that the BLR gas is virialized under the control 
of gravity of the central supermassive black hole, black hole 
mass Mbh can be estimated using the so-called linewidth- 
luminosity-mass scaling relation (e.g., Kaspi et al. 2000). 
We calculate the black hole masses using the formalism pre- 
sented in Vestergaard & Peterson (2006, their equation 5) 
with H/3 FWHM measured from the best-fit model of H/3 
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broad-line profiles. This formalism was calibrated with re- 
verberation mapping-based masses and used the R-L re- 
lation by Bentz et al. (2006) corrected for host galaxy 
starlight contamination. The Eddington ratios (I/boi/^Edd) 
are calculated assuming a bolometric luminosity correction 
Lboi ~ 91/5100 as for normal QSOs (Kaspi et al. 2000, Elvis 
et al. 1994). No correlation between the Balmer decrement 
and 1/5100, I/boi/^/Edd and Mbh are found to be significant 
at the Pnuii < 0.01 level (see Table 3). To check whether the 
results are dependent on the exact formalisms for black hole 
mass estimation, we re-calculate black hole masses using var- 
ious formalisms, including Peterson et al. (2004) and Onken 
et al. (2004) and the R-L relation from Kaspi et al. (2005) 
and Bentz et al. (2006), respectively. For the line width, we 
also use the second moment of the line profile (crnne, often 
referred to as "line dispersion") measured from the model 
line profile. To avoid possible contamination to the contin- 
uum luminosity from jet and host-galaxy starlight, we also 
try to use the H/? luminosity to calculate the black hole 
masses. Yet no correlations are found, either. We also con- 
sider the possibility that these correlations may be present 
but differ for different black hole masses, and as such the 
trend of the correlations could be reduced when the sample 
spans a large range of black hole masses. We divide the sam- 
ple into several logAfan bins with a bin size of 0.4 dex for 
each bin. We find that, in none of the bins, the correlations 
of Ha''/H/3'' with I/5ioo or I/boi/-f/Edd are present. The rela- 
tions of Hq:''/H/3'' with 1/5100 and Z/boi/-f/Edd are displayed 
in Fig. 4a and 4b, respectively, where the 122 objects in 
the Mbh bin of lO'' *^ - 10**'^ Mq are denoted as solid black 
squares. 

The luminosity I/5ioo ranges from 2 x 10"*^ to 2 x 10"'^ 
ergs with a median of 2 x lO^** ergs s^^. We compile 
two subsamples with L5100 > 5 X ergs s~^ ('high-L' 
sample, 45 objects) and Lsioo < 1 x 10*^ ergs s~^ ('low-L' 
sample, 49 objects), respectively. The mean and standard 
deviation of log(^^) in the 'high-L' subsample are 0.496 
and 0.046, respectively; and 0.479 and 0.049 respectively in 
the 'low-L' subsample. The mean values of the two samples 
are consistent with each other, with a chance probability of 
9 per cent according to the Student's f-test. The Eddington 
ratios range from 0.005 to 1.2 with a median of 0.2. We also 
compile two subsamples of objects with Lboi/I/Edd > 0.5 
('high-state', 57 objects) and I/boi/^/Edd < 0.03 ('low-state', 
31 objects), respectively, which have the mean and standard 
deviation of log(§|^) of 0.488 and 0.043, and 0.499 and 
0.058, respectively. Again, the two subsamples have mutu- 
ally consistent mean values, at a chance probability of 35 
per cent given by Student's t-test. 



3.2.2 Radio loudness 

We cross-correlate the sample with the Faint Images of the 
Radio Sky at Twenty cm (FIRST; Becker et al. 1995) Survey 
following the procedure described in Lu et al. (2007). For un- 
resolved FIRST sources a matching radius of 3" is used; for 
resolved sources the FIRST images are visually inspected. 
In this way 68 matches are obtained. For these matches, we 
define and calculate the radio-to-optical fiux ratio (radio- 
loudness) following Ivezic et al (2002), Ri = log{^^2f^), 
where fi and /20cm are the fiux densities at i-band and 



20 cm, respectively. For the rest 336 objects that were cov- 
ered by FIRST, we calculate the upper limit of Ri by tak- 
ing the FIRST detection limit of 1 mjy as the flux lim- 
its at 20 cm. Fig. 4c shows the Hq''/H/3'' ratio versus -Ri. 
The correlation is weak at most, as indicated hy p — 0.205 
incorporating both detections and upper limits of Ri. We 
also calculate the k-corrected radio power at 20 cm as 
PiOcm = 47rD|,/i„i/(l + z)^"'""'', where the radio spectral 
index a,. [Fv oc !/""") is assumed to be —0.5 for all the ob- 
jects. Similarly, there seems at most a weak correlation, as 
indicated by p = 0.209. 

It is well known that the radio-loudness distribution ap- 
pears to be bimodajfl with ~90% being radio-quiet (RQ) and 
~10% radio-loud (RL) (Kellermann et al. 1989, Ivezic et al 
2002 and references therein). Hence for the two populations 
the central engines or physical processes related to accre- 
tion/jet may be different (e.g., Blandford et al. 1990, Boro- 
son 2002, Falcke et al. 1996). To further examine whether 
the Balmer decrements are systematically distinct for the 
two populations, we compile a RL and a RQ subsample fol- 
lowing Ivezic et al (2002): 19 sources with Ri > 1 axe re- 
garded as RL, while the rest 385 sources covered by FIRST 
are RQ. The mean and standard deviation of log(-^2j.) in 
the RL subsample are 0.528 and 0.057, respectively; and in 
the RQ subsample they are 0.483 and 0.046, respectively. 
The mean Balmer decrement appears to be slightly larger 
in the RL subsample by 0.05 dex than in the RQ one. This 
difference is significant according to Student's t-test, with a 
chance probability <C 0.01. 

3.2.3 Line Profile 

It has been known for a long time that the Balmer decrement 
is larger in the line core than in the line wings (e.g., Shuder 
1982, Crenshaw 1986, Stirpe 1991, Korista & Goad 2004 
and references therein). As a result the integrated Balmer 
decrement might depend on the line profile to some extent. 
To investigate this issue, we compute the skewness (the 3rd 
moment) and kurtosis (the 4th moment) based on the best- 
fit model of H/3 broad-line profiles, as well as FWHM and 
(^linc obtained in §3.2.1. The correlation between the inte- 
grated Balmer decrement and skewness is at most weak, if 
exist at all, as indicated by p = —0.203. No correlations 
are found between the integrated Balmer decrement and 
kurtosis, FWHM and (Tiinc, respectively. Considering that 
skewness and kurtosis are sensitive to the errors in the line 
wings caused by the substraction of the Fell, [O HI] and 
Hell emission, we also compute some empirical yet robust 
profile parameters based on the model profile. We define 
three indices to characterize the asymmetry (viz. skewness, 
Al), shape {SI) and kurtosis {KIl) as, Al = (C(|) - 
C(i))/FWHM, 57 = {FW{\) + f W(|))/(2 x FWHM), 
and KIl = FW{l)/FW{\), following De Robertis (1985), 
Boroson & Green (1992) and Marziani et al. (1996), respec- 
tively. Where C(^),C(|) is the centroid at j, | maximum, 
and FW{\),FW{j), the full width at i,| maximum, re- 
spectively. These three dimensionless parameters thus de- 

* However, the genuineness of this bimodality is still a matter of 
debate; see, e.g., Hooper et al. (1995), White et al. (2000) and 
Cirasuolo ct al. (2003). 
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fined are not afi'ected by the choice of the rest frame. Colhn 
et al. (2006) characterize the broad-hne profiles by the ratio 
of FWHM to aiine, based on which the broad- line profiles 
of AGN are separated into two categories: the first, having 
^^."^^ < 2.35, are narrower lines with relatively extended 

wings; the second, having ^^^^^^ ^ 2.35, are broader lines 
being relatively flat-topped. This linewidth ratio is similar 
to the above kurtosis index (KIl), and is therefore also com- 
puted and denoted as KI2. We perform correlation analysis 
of the integrated Balmer decrement with the above various 
parameters, and find no significant corrections (see Table 3). 
The relation between the integrated Balmer decrement and 
KIl is displayed in Fig. 4d. 

Among our blue AGN sample, there are 23 objects hav- 
ing Ha or H/3 lines of double or even multiple peaked pro- 
file^ (Shang et al., in preparation), with 5 being radio-loud. 
Such profiles are found in about 3 per cent of AGN (Strat- 
eva et al. 2003) and have a higher occurrence in LINERs 
(Eracleous 2004) and in RL sources (~20%, Eracleous & 
Halpern 2003). Using a sample selected from RL AGN, Era- 
cleous & Halpern (1994) found that double-peaked emitters 
have the large integrated Balmer decrements, like their par- 
ent population of RL AGN; later they further found that 
their integrated Balmer decrements are even larger than 
that of the latter (on average 5.23 versus 4.26, Eracleous 
& Halpern 2003). We investigate this issue using the above 
23 objects (DBF subsample). The mean and standard devi- 
ation of log(^^) in the DBF subsample are 0.515, 0.055, 
respectively, while these values are 0.484 and 0.046 in the 
rest of objects (non-DBF subsample). The mean integrated 
Balmer decrement appears to be slightly larger in the DBF 
subsample than in the non-DBF sample, which is signifi- 
cant with a chance probability of 0.002 (0.01) by Student's 
t-test assuming the two distributions to have the same (dif- 
ferent) variance. We cannot find any difference in integrated 
Balmer decrement between the DBF subsample and the RL 
subsample (with a f-test probability of 44%), possibly due 
to the small size of the two subsamples. 



3.2.4 Other SED properties 

The shape of the ultraviolet to X-ray ionizing continuum 
affects the extended partially ionized zone and thus the 
Balmer decrements (Kwan & Krolik 1981). Here we take the 
Qox values (the slope of a hypothetical power law between 
25OOA and 2 keV) for 268 matched sources from Anderson et 
al. (2007). These values were derived from broadband X-ray 
(ROSAT) and g-band (SDSS) fluxes by assuming an X-ray 
energy index Ox ~ 1.5 and an optical index ax = 1.5. No 
correlation between the Balmer decrement and «ox is found. 
We also perform correlation analysis between the Balmer 
decrement and the optical-near-ultraviolet continuum slope 
(oa.i); and flnd no significant correlation (see Table 3). The 
relations of the Ha''/H/3'' ratio with aox and Oa,! are dis- 
played in Fig. 4e and 4f, respectively. 



® Such profiles are defined as having the number of peak and 
"pseudo-peak" greater than 2; a "pseudo-peak" is defined as the 
point where the 2nd derivative is minimal and negative. See Shang 
et al. (in preparation) for details. 



4 DISCUSSION 

4.1 Representativeness of Unreddened AGN 

Internal Dust extinction in objects of our blue AGN sample 
should be negligible. In the sample, qa,i varies in the range 
from 1.5 to 2.7 (see Table 1) with a standard deviation of 0.2. 
Such a scatter is consistent with previous results that there 
is a large intrinsic dispersion of the continuum slope (Elvis 
et al. 1994, Rowan-Robinson 1995, Natali et al. 1998, Kuhn 
et al. 2001). To check how much the sample suffers from ex- 
tragalactic dust extinction, we cross-match these AGN with 
the Galaxy Evolution Explorer (GALEX; Morrissey et al. 
2005) General Data Release 3 with a matching radius of 
2. "6 in a way similar to Trammell et al. (2007). We obtain 
253 objects with both reliable far-ultraviolet (/) and near- 
ultraviolet (n) magnitudes, while 185 of the rest have not 
been covered by GALEX yet. For the 253 objects we com- 
pute the relative f—n colors defined as the f — n colors minus 
the corresponding median color of quasars at the same red- 
shift. The digital curve of median color versus redshift is 
kindly provided to us by G. Trammell, as presented in Fig. 
13 of TrammeU et al. (2007). We find the distribution of the 
relative f — n colors can be well described by a Gaussian 
with a mean of zero and a standard deviation of 0.3, a dis- 
persion caused by GALEX photometric uncertainties only 
(TrammeU et al. 2007, see their Fig. 14). We also find that 
neither the f — n color nor the relative f — n color correlates 
with the Balmer decrement. Hence we believe that there is 
little dust extragalactic extinction in our sample. 

We notice that there exists a potential bias against se- 
lecting AGN with intrinsic red slopes by our criterion. To 
test the representativeness of our sample, we perform Spear- 
man correlation analysis of the continuum slope with all the 
various, above-mentioned parameters in the same way as in 
§3.2. The results are listed in Table 3. It turns out that no 
correlations of aA,i are present with the Eddington ratios, 
black hole mass, radio-loudness, radio power, linewidth and 
various other line-profile parameters. But there is a positive 
correlation between the continuum slope and the nuclear lu- 
minosity L5100 as indicated by p = 0.457 and P„uii = 10~® 
(see Table 3)0 Thus our sample is biased toward relatively 
high luminosity. This is consistent with the fact that the 
sample covers a high luminosity range as L5100 Si 2 x 10'''' 
and Hq luminosity > 5 x 10'*'^ ergs s^'^. To sum up, our sam- 



" Such a slope— luminosity correlation has also been reported in 
the literature (e.g., O'Brien et al. 1988, Cheng et al. 1991, Francis 
1993, Carballo et al. 1999, Kuhn et al. 2001). Previously, how- 
ever, it was not clear whether the slope— luminosity correlation is 
intrinsic or just arises from the correlation between the slope and 
redshift in a magnitude-limited sample. Some authors claimed an 
intrinsic dependence of the continuum slope with redshift rather 
than with luminosity (O'Brien et al. 1988, Cheng et al. 1991, 
Francis 1993). Our sample is much larger than the previous ones 
and has a limited redshift range z < 0.35. In our sample no sig- 
nificant correlation between the continuum slope and redshift is 
found; yet the significant correlation between slope and L5100 still 
exists by partial correlation analysis under the control of redshift, 
with a coefficient of 0.36 and a chance probability <C 10~^. As 
a by-product of this study, we suggest that the slope-luminosity 
correlation is intrinsic. 
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pie is free from biases regarding the above-mentioned AGN 
properties except luminosity. 

4.2 On the Intrinsic HaVH/J' Distribution and Its 
Applications 

As discussed in §1, historically the broad- line Balmer decre- 
ments in Seyfert Is and QSOs have been generally thought 
to be considerably steeper than the Case B value and to 
have a large intrinsic dispersion. Such a belief may have 
originated from results of early studies with small samples, 
and was (mis-)reinforced by observations of Seyfert 1.8/1.9 
and broad-line radio galaxies, and further assured by its con- 
sistence with photoionization modeling mainly of individual 
clouds in various conditions (see observational and theoreti- 
cal references in §1). In a sample of 36 Seyfert Is (including 
Seyfert 1.2 and 1.5) studied by Osterbrock (1977), the mean 
Hq/H/3 ratio is 3.6 with a range from 2.6 to 5.9 (uncorrected 
Galactic reddening). Similar results were found in a sample 
of 24 Seyfert 1/QSOs by Neugebauer et al. (1979), giving 
the mean Ha/H/3 ratio 3.6 with a range from 2.2 to 4.9 
(uncorrected Galactic reddening; Ha measurements include 
[Nil] flux). In addition to relatively large errors inherent in 
the early measurements and the neglect of Galactic redden- 
ing correction, we believe that internal reddening also plays 
a non-negligible role in causing the steeper Balmer decre- 
ments in previous results. Similar situation was found in 
PG quasars (Rowan- Robinson 1995, Baskin & Laor 2004). 
For instance, I ZW 1 has an observed Ha/H/3 ratio of 4.86 as 
given in Table 1 of Osterbrock (1977), but actually its broad 
lines suffer internal reddening of E{B — V) ^ 0.1 (Rudy et 
al. 2000) as well as Galactic reddening of E{B - V) ~ 0.1 
(Schlegel et al. 1998; Stark et al. 1992). Interestingly, for the 
9 Seyfert Is listed in Table 2 of Wu et al. (1983), the mean 
Ha''/H/3'' ratio is 3.1 after individual correction for Galactic 
reddening (Goodrich et al. 1990), which is almost the same 
as ours found here. Recently, in a sample of broad-line AGN 
with low starlight contamination, Greene & Ho (2005) found 
that the mean value of the total (narrow + broad) Ha/H/3 
ratio is 3.5; the somewhat larger value compared to ours 
should also be due to internal reddening. In a large sample 
of about 2000 narrow-line Seyfert Is (NLSls), Zhou et al. 
(2006) found the mean Ha''/H/3'' ratio to be 3.0, close to the 
value derived here for blue AGN. In summary, results ob- 
tained from early small samples and from recent much large 
samples all point to that the mean Ha''/H/3' ratio of Seyfert 
1/QSOs is only slightly larger than the Case B value. 

A surprising feature emerged from the present study is 
the actual dispersion of the intrinsic Ha''/H/3'' ratio being 
rather small (0.03 dex in a log-normal distribution)! This is 
contrary to the prevalent belief that there is a "considerable 
range of intrinsic (Balmer) line ratios" (Wu et al. 1980), 
which was based on early samples that were in fact too 
small to make a general conclusion. Our conclusion should 
hold for luminous Seyfert Is and QSOs in general, since by 
selection our sample is only biased against AGN with in- 
trinsic red continua, but the Balmer decrement does not 
correlate with the continuum slope as found in §3. We fur- 
ther suggest that our result is also likely to hold for AGN 
of low luminosity. Firstly, there is no correlation between 
the Balmer decrement and luminosity in our sample, and 
the 'high-L' and 'low-L' subsamples have the indistinguish- 



able mean Balmer decrements. Secondly, in the large sample 
of about 2000 NLSls that are generally at high accretion 
states, the Balmer decrements still cluster tightly around 
~3 even when the nuclear Lsioo goes down to 10*^ ergs s~^ 
(Zhou et al. 2006). Moreover, as shown in §3, the Balmer 
decrement does not correlate with the Eddington ratio. Re- 
cently, contrary to our findings here. La Mura et al. (2007) 
found a weak correlation between the Ha''/H/3'' ratio and 
the Lboi/^Edd- However, as suggested by those authors, this 
might result from the inclusion of reddened objects in their 
sample. They also found a weak correlation between the flux 
ratios of broad-line Balmer series (up to H(5) and the line 
width; but, in the case of Ha''/H/3'' ratio, the correlation is 
not significant as indicated by Pnuii = 3.34 x 10~^. 

It should be noted that in some variable objects an 
anti-correlation between the continuum flux (i.e., accretion 
state) and the Ha''/H/3*' ratio has been reported during flux 
variability over a time scale of months (e.g., the prototy- 
pal NGC 5548; Wamsteker et al. 1990, Dietrich et al. 1993, 
Shapovalova et al. 2004). Such a behavior has been well 
predicted and explained by theoretical models (e.g., Net- 
zer 1975; Korista & Goad 2004). Thus the lack of such a 
correlation in the AGN ensemble seems to be somewhat un- 
expected. We guess that, instead of instantaneously respond- 
ing to the continuum variation, the BLR clouds distribution 
(as a function of mass and/or luminosity) is perhaps deter- 
mined/adjusted by the long-term average of the accretion 
state; for an AGN ensemble, this long-term adjustment may 
make the average Ha''/H/3'' ratio in a narrow range. In ad- 
dition, according to Table 8 of Shapovalova et al. (2004) 
and Table 2 of Dumont et al. (1998), the Balmer decrement 
of NGC 5548 varies between 3.0 and 4.3, which is not ex- 
treme and is actually within the range for the blue AGN 
ensemble as reported here. Due to the small variability am- 
plitude, any correlation, even if exists in individual objects, 
would be smeared out in the ensemble. In a few cases, such 
kind of variability can be explained by variation in extinc- 
tion (Goodrich 1989, Tran et al. 1992, Goodrich 1995; see a 
discussion about the time scale of this kind of variability in 
Wang et al. 2007). 

From theoretical perspective, however, a tight Ha''/H/3'' 
distribution around 3.1 seems to be quite surprising — all the 
BLR photoionization computations predict a rather large 
range in the Balmer decrement for individual clouds in plau- 
sible BLR conditions (see references in §1). It seems unlikely 
that such a discrepancy is mainly caused by the incapability 
of photoionization modeling of the BLR clouds, because all 
the model computations over the past 40 years gave the simi- 
lar trend of large HaVH/3'' range (e.g., Netzer 1995, Dumont 
et al. 1998, Korista & Goad 2004), although such computa- 
tions for a single cloud are still uncertain (see Netzer 1995 
and references therein). Another possibility, in fact an old 
proposal, is that the BLR emitting gas of all AGN have been 
"fine-tuned" to a certain ionization parameter by some phys- 
ical processes (e.g., the "hot- warm" model of Krolik et al. 
1981). However, the existence of such a fine-tuning, among 
clouds in likely jumbled environments as in AGN BLRs, is 
questioned by others and deemed to be unnatural and un- 
likely (e.g., Mathews & Ferland 1987; Baldwin et al. 1995). 
Here, we suggest a plausible explanation to this seemingly 
discrepancy invoking the "locally optimally-emitting cloud" 
(LOG) model, as proposed by Baldwin et al. (1995). The 
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essential idea of the LOG model is that each line arises pre- 
dominantly from clouds only in a narrow range of density 
and distance from the continuum source, due to natural se- 
lection effects largely introduced by the atomic physics. In 
this scenario, the similar values of the Balmer decrements, 
just like other surprising similarities in emission-line spec- 
tra of Seyfert 1/QSOs (Davidson & Netzer 1979, Baldwin et 
al. 1995), appear to be a natural consequence. According to 
the recent calculation by Korista & Goad (2004, their Fig- 
ure 5) using the spectral synthesis code CLOUDY (version 
90.04), the HaVH/3'' flux ratio varies from roughly 17 to ap- 
proximately 1 across the parameter plane of gas density and 
ionizing flux characteristic of BLR clouds; the Ha''/H/3'' ra- 
tio, when integrated over the full BLR of their LOG model, 
varies from 3.7 (high state) through 4.9 (low state). Using 
GLOUDY version 07.02 (last described by Ferland et al. 
1998) with improved collisional rates for excited states of hy- 
drogen atoms, the Ha''/H/3'' ratio contours now have values 
that are typically ~0.1 dex sma//er than those presented in 
Figure 5 of Korista & Goad (2004) (Korista, private commu- 
nication) . Thus there exists a large span of cloud parameters 
within the density-flux plane, over which most of the broad 
emission lines are emitted, for which the Hq''/H/3'' ratios lie 
between 2.5 and 4. In fact LOG integrations over the range 
of cloud parameters adopted by Korista & Goad (2004) now 
predicts an Ha^'/H/S'' ratio that is consistent with our mea- 
surements (Korista, private communication). We therefore 
suggest that the LOG model, together with improvement 
on the photoionization modeling, can give a natural expla- 
nation to the seemingly discrepancy between the theory and 
the measurements. 

As is found here, radio-loud objects and objects having 
double-peaked emission-line profiles among the blue AGN 
ensemble have slightly larger Balmer decrements than the 
rest on average. In these kinds of AGN the physical condi- 
tions related to accretion/jet may be different (Blandford 
et al. 1990, Falcke et al. 1996, Eracleous & Halpern 2003), 
e.g., with additional X-ray ionizing radiation from the jet 
base. Thus collisional excitation, self-absorption, and other 
line-transfer effects may play roles to enlarge the Balmer 
decrements, as discussed by some authors (e.g., Osterbrock 
et al. 1976, Grenshaw et al. 1988, Netzer et al. 1995). Given 
the small fractions of these two classes (~ 10 and 3 per cent, 
respectively) among the AGN population, the distribution 
of the Ha''/H/3' ratios show little changes when the RL and 
double-peaked objects are excluded. In another word, for the 
bulk of the AGN population, such effects are insignificant. 

Our finding out of this study may have an interest- 
ing implication: the precise distribution of the intrinsic 
Hq^/H/S*' ratio, that is insensitive to any nuclear properties 
known so far, renders it a useful tool with which the BLR 
extinction can be derived, at least in a statistical manner. 
Specifically, for a sample concerned, we can derive its in- 
ternal E{B — V) distribution by de-convolving the observed 
HaVH/3'' distribution with the intrinsic HqVH/?'' distribu- 
tion as found here. This approach has now been applied to a 
large sample of broad-line AGN culled from the SDSS DR4 
to derive the internal E{B—V) distribution of the AGN BLR 
in the local universe, and hence to obtain the fraction of ob- 
scured AGN of various intrinsic luminosity and of various 
degree of extinction (Zhang et al., in preparation). At the 
zeroth-order approximation, Dong et al. (2005) have used 



the Balmer decrement to derive the BLR reddening for a 
sample of luminous Seyfert 1.8/1.9 galaxies culled from the 
SDSS EDR and derived the fraction of partially obscured 
quasars in the local universe. Follow-up XMM-Newton ob- 
servations confirm that these sources have large absorption 
column densities in the X-ray (Zhou et al., in preparation). 



5 SUMMARY 

We have investigated the broad-line Balmer decrements for 
a large, homogeneous sample of 446 blue AGN, of Seyfert 1 
galaxies and QSOs. They are selected from the Sloan Dig- 
ital Sky Survey Fourth Data Release according to the cri- 
teria of redshift z < 0.35, the median spectral signal-to- 
noise per pixel ^ 10, and the continuum slopes qa ^ 1.5 
[fx ~ A~°*) that are fitted in the rest-wavelength range 
of 4000-5600 A. With the blue criterion of the continuum 
slope, dust extinction in the sample objects is expected to 
be negligible, which is also confirmed by their relative col- 
ors in the ultraviolet. The sample is fairly representative of 
normal Seyfert 1/QSOs (at least the luminous objects with 
Ha luminosity greater than lO**^ ergs s~^), in light of the 
fact that the optical-near-ultraviolet continuum slope does 
not correlate with the Balmer decrement, nor with other 
AGN properties except nuclear luminosity. We find that (i) 
The distribution of the intrinsic broad- line Ha/II/3 ratios 
can be well described by log-normal, with a peak at Ha/II/3 
=3.06 and an dispersion of likely 0.03 dex only, (ii) there 
are no significant corrections between the Balmer decre- 
ment and the nuclear properties such as luminosity, accre- 
tion rate, continuum slope and aox', (iii) on average, the 
Balmer decrements are slightly larger in radio-loud objects 
(3.37) and objects having double-peaked emission-line pro- 
files (3.27). Therefore we suggest that the broad-line Ha/H/3 
ratios can be used as a good indicator of dust extinction of 
the AGN broad-line regions, at least in a statistical manner. 
This result is especially true for radio-quiet AGN with reg- 
ular emission-line profiles that constitute the vast majority 
of the AGN population. Such an application has significant 
implications for deriving the distribution of internal dust 
extinction in the BLR of AGN, and hence the obscuration 
fraction of AGN in the universe. 
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Figure 1. Representative examples of the SDSS spectra and simultaneous fitting of the continuum, Fell and other emission lines from 
4030A to 75OOA. In each panel, we plot the observed spectrum (black), the sum of the best-fit components (red), the continuum modeled 
as a broken power-law with a break at 5600A (blue), and the Fell emission (green). The radio-loudness {Ri = log( ^^Pg"* )), the continuum 

slopes blueward {ax 1) and redward {ax 2) of 5600A are indicated. Note that the model Balmer broad-line profile in some objects differ 
more or less from the observed one due to that the Balmer lines are assumed to have the same profile in the fit; for such objects we refit 
the line profiles to obtain much accurate line parameters (see the text in §2.4 and Fig. 2). 
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Figure 2. Results of our line-profile fitting procedure applied to the H/3 region (left panels) and the Ha region (right panels) for the 3 
objects demonstrated in Fig. 1. We plot the original data (black), the sum of all the best-fit components (red), the fitted narrow lines 
(green), the fitted broad H/3 and Ho (blue), and the residuals of the fit (bottom, offset downward for clarity). The broad-line Balmer 
decrements are indicated. 
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Figure 3. Histogram of log( ) for the 446 blue AGN. Also displayed are the fit (solid line) with a Gaussian function yielding a mean 

of 0.486±0.002 and a standard deviation of 0.046±0.002. The intrinsic log(-|^) distribution is over-plotted (dash line) with a standard 
deviation of 0.03, that is estimated by de-convolving the observed distribution with the dispersion caused by measurement uncertainty. 
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Figure 4. Plots of HqVH/?'' versus ALa(5100A), Lbol/^'Eddi radio loudness Ri, kurtosis index KIl as defined by Marziani ct al. (1996), 
c«ox, and the optical— near-ultraviolet continuum slope a\. In panel (a) and (b), objects in the Mbh bin of 10^'* — 10*^'^ Mq are denoted 
as solid black squares. In panel (c), the 404 objects covered by FIRST are plotted, among which 336 have only upper limits (grey circles 
with arrows). In panel (e), the 268 sources having ROSAT matches are plotted, where the Oox values are taken from Anderson et al. 
(2007). Note that although the moan Balmer decrements are slightly larger in the radio-loud subsample and the subsamplc of objects 
with double-peaked line-profiles, there are no significant correlations of the Balmer decrement with Ri or KIl in the blue AGN ensemble. 
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Table 1. Properties of Blue AGN 



Object 






J^sioo 


F(Hf)") » 


F(H/3'') " 


F(Ha") ° 


F(H<:<'>) " 


°A,1 


°A,2 


/ 








J20CTT1 


(1) 




(2) 


(3) 


(4) 


(6) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


000710.01 + 005329.0 





31620 


3.85 


27.8 


2778.8 


146.0 


9494.1 


1.83 


1.83 


18.24 


17.66 


17.39 


0.032 


1.44 


000834.72 + 003156.1 





26303 


2.71 


191.4 


1002.1 


574.2 


2590.6 


2.02 


1.61 


18.46 


17.95 


17.48 


0.040 


<1 


000943.14-090839.1 





20958 


2.09 


6.7 


2932.7 


41.7 


7391.5 


1.57 


0.92 






16.92 


0.038 


<1 


001224.02-102226.2 





22822 


2.69 


52.3 


3510.8 


147.0 


11071.3 


1.81 


1.10 


17.79 


17.83 


16.61 


0.037 


3.73 


001247.93-084700.4 





22006 


1.95 


0.0 


2676.2 


61.8 


8422.0 


1.70 


1.39 






16.92 


0.032 


<1 


002840.69-102145.0 





32173 


1.53 


4.2 


493.6 


27.0 


1608.8 


1.58 


1.35 






18.70 


0.037 


<1 


004319.74+005115.3 





30807 


3.08 


37.7 


2149.9 


176.6 


7098.0 


1.72 


1.72 


18.21 


18.08 


17.88 


0.020 


1.60 



Note. — Col. 1, object name in J2000.0. Col. 2, redshift given by the SDSS spectroscopic pipeline. Col. 3, monochromatic 
luminosity \L\ at 5100 A, In units of ergs s~^. Col. 4, H/3 narrow component flux; its typical error is 6%. Col. 5, H/3 broad 
component flux; its typical error is 8%. Col. 6, Ha narrow component flux; its typical error is 5%. Col. 7, Ha broad component 
flux; its typical error is 5%. Col. 8-9, continuum slopes blueward and redward of 5600A, respectively (/a = A~"-^); their 
typical error is 8%. Col. 10-11, GALEX calibrated magnitudes (AB) in the FUV (/) and NUV (n) bands, respectively, 
uncorrected for Galactic extinction; a "—999" is given for sources that are covered yet not detected by GALEX, and a blank 
is given for 185 sources that are not covered by GALEX. Col. 12, the SDSS i-band magnitude (AB), uncorrected for Galactic 
extinction. Col. 13, the Galactic color excess derived from Schlegel et al. (1998). Col. 14, the integrated flux density at 20cm 
detected by FIRST, in unit of mjy; the detection limit of 1 mjy is adopted as the upper limit for 336 objects that are covered 
yet not detected by FIRST, and a blank is given for 42 objects that are not covered by FIRST. Table 1 is now available in 
its entirety at \http: //staff. ustc. edu.cn/~xbdong/Data_Release/blue AGN_DR4/\ it will be available via the link 
to the machine-readable table on the MNRAS website. 

=1 Line flux in units of 10~^^ erg s~^ cm~^. 



Table 2. Ha and H/3 Broad-Line Profile Measurements 







broad H/3 








broad Ha 






Object 


FWHM 


0"liiio 


Al 


SI 


KIl 


FWHM 




Al 


SI 


KIl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


000710.01+005329.0 


9165.5 


4798.8 


-0.16 


0.93 


0.38 


9165.5 


4798.8 


-0.16 


0.93 


0.38 


000834.72+003156.1 


1953.0 


1415.8 


0.07 


1.08 


0.40 


1953.0 


1415.8 


0.07 


1.08 


0.40 


000943.14-090839.1 


5288.3 


3406.5 


-0.04 


1.05 


0.47 


4225.2 


2236.0 


-0.01 


1.05 


0.53 


001224.02-102226.2 


4294.2 


4175.2 


-0.64 


1.40 


0.18 


4710.4 


3646.0 


-0.35 


1.15 


0.26 


001247.93-084700.4 


3007.4 


2488.0 


-0.10 


1.13 


0.35 


3007.7 


2488.0 


-0.10 


1.13 


0.35 


002840.69-102145.0 


1371.8 


1525.8 


-0.05 


1.13 


0.36 


1371.8 


1525.8 


-0.05 


1.13 


0.36 


004319.74+005115.3 


12630.6 


5827.4 


-0.02 


1.01 


0.71 


12630.6 


5827.4 


-0.02 


1.01 


0.71 



Note. — Col. 1, object name in J2000.0. FWHM and line dispersion (aiino) in unit of km s~^. 
Asymmetry index (Al) defined as in De Robertis (1985); shape index (SI) as in Boroson 
& Green (1992); kurtosis index (KIl) in Marziani et al. (1996); see §3.2.3. All parame- 
ters are derived from the model broad lines. Table 2 is now available in its entirety at 
http://staff.ustc.edu.cn/~xbdong/Data_Release/blueAGN_DR4/; it will be available 
via the link to the machine-readable table on the MNRAS website. 
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Table 3. Summary of Spearman Rank Correlation Tests ^ 





net / rip 




Z/5XQQ 


0.115 (0.015) 


0.457 (1.0 X 10-*^) 


D OI / H; Q Q 


-0.051 (0.281) 


0.073 (0.123) 




0.093 (0.049) 


0.187 (6.8 X 10^5) 


Ri " 


0.205 (< 10""') 


0.044 (0.381) 




0.209 (< 10"") 


0.144 (0.004) 


FWHM 


0.078 (0.098) 


0.075 (0.116) 




0.099 (0.037) 


0.121 (0.011) 


skcwncss 


-0.203 (1.5 X 10-5) 


-0.077 (0.103) 


kurtosis 


-0.047 (0.324) 


-0.055 (0.242) 


AI d 


-0.020 (0.673) 


0.024 (0.619) 


SI 


-0.002 (0.974) 


0.129 (0.006) 


KIl <i 


-0.151 (0.001) 


-0.070 (0.138) 


Kn d 


0.047 (0.320) 


-0.032 (0.502) 


"ox ° 


-0.089 (0.147) 


0.138 (0.024) 




-0.109 (0.022) 





For each entry, we list the Spearman rank correla- 
tion statistic (p) and the probability of the null hypoth- 
esis (Pnuii) in parenthesis. If no censored data (upper 
limits) are present, p is equal to the Spearman's rank 
correlation coefficient Vs- 

^ The black hole masses are calculated using the 
formalism presented in Vestergaard & Peterson (2006, 
their equation 5) with H/3'' FWHM listed in Table 2; 
Eddington ratios (Lboi/I'Edd) are calculated assuming 
that the bolometric luminosity Li,oi ~ 9L5100 used for 
normal QSOs (Kaspi et al. 2000, Elvis et al. 1994). 

Using 404 sources covered by FIRST, of which 336 
are upper limits. 

Computed based on the model H/3 broad-line pro- 
files. 

" Using 268 sources matched with ROSAT; Qox values 
are brought from Anderson et al. (2007). 
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